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Abstract. The interaction of a two- level atom with a single- mode quantized field is 
one of the simplest models in quantum optics. Under the rotating wave approximation, 
it is known as the Jaynes-Cummings model and without it as the Rabi model. 
Real-world realizations of the Jaynes-Cummings model include cavity, ion trap and 
circuit quantum electrodynamics. The Rabi model can be realized in circuit quantum 
electrodynamics. As soon as nonlinear couplings are introduced, feasible experimental 
realizations in quantum systems are drastically reduced. We propose a set of two 
photonic lattices that classically simulates the interaction of a single two-level system 
with a quantized field under field nonlinearities and nonlinear couplings as long as the 
quantum optics model conserves parity. We describe how to reconstruct the mean 
value of quantum optics measurements, such as photon number and atomic energy 
excitation, from the intensity and from the field, such as von Neumann entropy and 
fidelity, at the output of the photonic lattices. We discuss how typical initial states 
involving coherent or displaced Fock fields can be engineered from recently discussed 
Glauber-Foek lattices. As an example, the Buek-Sukumar model, where the coupling 
depends on the intensity of the field, is classically simulated for separable and entangled 
initial states. 



PACS numbers: 05.60. Gg, 42.50.Ex, 42.79.Gn,42.82.Et 

OCIS numbers: (350.5500) Propagation; (230.4555) Coupled Resonators; (230.5298) 
Photonic Crystals; (270.0270) Quantum Optics; (270.5580) Quantum Electrodynamics; 
(310.2785) Guided Wave Applications. 



1. Introduction 



Optical analogy has been of great importance to the development of quantum mechanics; 
cf. [l]. Quantum information processing has been one of the fields where authors have 
extensively studied classical analogies by use of linear optics elements involving classical 
light beams [2] or single-photons [3]. In recent years, the advent of technology that 
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allows manufacturing arrays of optical waveguides has detonated the classical simulation 
of quantum and relativistic systems [IHIH]. 

We are interested in the classical simulation of a single two-level atom interacting 
with a quantized field, the basic building block in quantum optics, under nonlinear 
processes in the field and couplings as the linear interaction has already been classically 
simulated [191 120]. The quantum model of interaction between an ensemble of two- 
level atoms and a radiation field under a series of approximations was introduced 
by Dicke |21]. The exact solution for the single atom case was found by Jaynes and 
Cummings [22] and the Jaynes- Cummings (JC) model became a work horse of quantum 
optics as well as the Dicke model. Among other interesting phenomena, the JC model 
shows a collapse and revival of the mean atomic excitation energy when interacting 
with coherent radiation [23], but it has not been possible to describe it in closed form. 
Buck and Sukumar (BS) presented a JC model with nonlinear coupling where the mean 
atomic excitation energy can be evaluated in exact closed form and exhibits a periodic 
collapse and revival similar, for small times, to that seen in the JC model [21] . The most 
general JC model involving nonlinear couplings, up to our knowledge, was presented and 
solved by Kochetov [25] . 

H = iOfa^a + ujqS, + A \a^''a^ f{a^a)S+ + f{a^a)d^^a''sJ\ . (1) 

The Kochetov model describes the interaction between a single atomic system, with 
a finite number of equidistant levels, and a single-mode radiation field. The atom 
is described by the generators Si with i = z,+,— that obey an su{2) algebra. 



±S± and 



= 2Sz, and the inter- level energy is given by ujq. The 

field is described by the creation (annihilation) operators, a^{a), and the frequency 
Uf. The parameter A is a coupling constant and the function f{a^a) is a real well- 
behaved function of the number operator. Peculiar phenomena has been found among 
the years with specific realizations of the Kochetov model; e.g. for just two levels [261432] . 
Some other nonlinear models studied along the years for the two-level case involve: 
driving [33], field nonlinearities [Mti36] . or consider the coupling without the rotating 
wave approximation (RWA) [37]. The great majority of these models does not have a 
feasible experimental realization with current quantum technology. 

In this contribution, we are interested in a specific class of the Kochetov model that 
we extend to consider field nonlinearities and the counter rotating terms neglected by 
the RWA. We present our nonlinear Rabi model and its characteristics in the following 
section; in short, it conserves parity and this allows us to split the Hilbert space in 
two parity basis. Then, we show that these two parity basis allow us to describe 
the dynamics of our quantum model with classical light propagating in a set of two 
photonic lattices. We show how to construct typical initial states; e.g. an initial state 
involving just a photon number state and the ground or excited state corresponds to 
light impinging just one waveguide and those involving coherent states can be prepared 
through Glauber- Fock photonic lattices [381 - HO] . We also discuss the most common 
quantum optics measurements and show how they are related to measurements of the 
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intensity at the end of the classical simulators, with just an exception that requires 
knowledge of the amplitude and phase at each output. In order to provide a working 
example, we present the BS model in Section 4 and close with a brief conclusion. 



2. The quantum optics model 

We are interested in a general nonlinear Rabi model, 

H = h[n) + —cTz + 9~ 0- — 7^^+ H 7^'^^ 



n 



+ 9.(aM*_ + Mat,,), (2) 



'n yn 

with well-behaved real functions h{n) and f{n) in terms of the number operator, 
n = a^a, and where the operators with i = z,+,— are Pauli matrices. We have 
split the nonlinear Rabi coupling f{h)/^/fi-\- f{n)/^/fi a^^ + a^) into Jaynes- 
Cummings coupling described by the coupling parameter g_ and counter-rotating terms 
described by the parameter g+ for reasons that will become apparent in the next section. 
As mentioned before, this general Hamiltonian is not physical realizable with current 
experimental setups. For example, in order to implement ([2]) in cavity- or circuit-QED 
the function h{n) must be linear to account for the free-field energy or quadratic at most 
to describe the effect of a Kerr medium, h{n) = ujfi + nn^, and the coupling function 
f{n)/y/n must be a constant. Spin chains [H], ion traps |12] - H^ . and atoms in optical 
lattices |15], just to mention a few examples, may provide a way to realize some of the 
Hamiltonians covered by this class as far as one could attain precise control of the size, 
self-energies and couplings in the spin chain analog, drive an ion with multiple lasers in 
the second, and engineer adequate optical lattices in the latter. 

The Hamiltonian in (|2]) has constant of motion: parity, which can be defined as 



that is 
bases. 



n = -{-If a,; (3) 
iy,n =0. Conservation of parity allows us to define two orthogonal parity 

|+,j)= (fity|0,^) = {|0,(7),|l,e), 12,(7),...}, (4) 

|-,j)= (5ty|o,e) = {|0,e), 11,(7), |2,e),...}, (5) 

with = (l/y/hTl)a)a^ jlHHlH]. These parity bases split the Hilbert space in two 
orthogonal parity subspaces which are the foundation for the classical simulation of our 
nonlinear Rabi model. 

In addition, if the counter-rotating terms coupling parameter is null, = 0, then 
the model reduces to a nonlinear JC model that also conserves the total number of 
excitations defined as = n + cFz/'^'-, i-e. H,N =0. It is simple to obtain the 
time evolution for this case by using a method that makes use of Susskind-Glogower 
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operators [H], V = J2T=o + 1|; such that (j2]) with g+ = becomes 

rr^ , V 0\( hin + l)+Uo/2 g_f{n) \ ( ^ . 

' 1 ; 1^ g^fin) hin + l)-uo/2 )[ 1 ' ' ^ ^ 

\> \ / T{h) + n{h) T{n)-n{h) \ f \+{n) 
Q I )\ 2g_f{n) 2g^f{n) J [ A_(n) 

T{h) + n{h) r(n) - y Y o\ 

2g^f{n) 2g.f{n) J [ 1 J' ^' ^ 

with the elements of the similarity transformation given by 

r{h) =h{h-l)-h{n)+uJo, (8) 



n{n) = ^T\h)+4glPin), (9) 
and the dispersion relation, 

XJn) = '''''-"+^^'">=^"''''. (10) 
Photon transport is then given by the time evolution operator, 

~ ~ 1 )[ 2g^f{h) 2g^f{n) ) [ e-^-(")* 

1^ 2g_f{h) 2g_f{n) J [ 1 J' ^ ' 

In other words, by using the analogy between transport of single-photon states and 
propagation of classical field, it is very simple to calculate the propagation through the 
equivalent nonlinear JC photonic lattice via quantum optics methods. 

3. Classical simulation in arrays of coupled waveguides 

The dynamics of a Hamiltonian system is given by the Schrodinger equation 

tdtim) = H\m), (12) 

where the notation dx stands for partial derivative with respect to x. As our model 
Hamiltonian ([2]) conserves parity, we can split the evolution in even and odd parts, 
IV') = \'4'^~^K't)) + l^*'"^^))' ^iid propose the solutions 

oo 

\^^^\t))=J2^i^\t)\±,j). (13) 

3=0 

This parity decomposition leads to one coupled differential set for each of the parity 
subspaces: 

td,£t^ = ^iW(0)£:f ) + g±f {l)£i^\ (14) 
td4tli = d^^H^k + l)£^tli + 9±f{2k + l)£^t^ + g^f{2k + 2)4^j2, k > 0, 
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(15) 

^^tS^t'' = d^^\2k)£^i^ + g^f{2k)£^il, + g^f{2k + l)£^£^,, k > 1, 

(16) 

with 

d^^\j) = Hj) t (-l)^y. (17) 

Each of these coupled differential sets is equivalent to that describing an array of 
waveguides coupled to their nearest neighbor up to a global phase and replacing t ^ z. 
In this analogy, the refraction index of the jth waveguide is modulated by d'^^\j) and the 
separation distance between the jth and the preceding waveguide is proportional to /(j). 
We are interested in these one- dimensional photonic crystals due to the great control that 
can be attained in modulating both the refractive index of individual waveguides and the 
separation between them by femtosecond laser waveguide writing on fused silica [50]. In 
recent experiments, the precise control over the refractive index and coupling parameters 
of photonic lattices has allowed the classical simulation of the Rabi model for coupling 
parameters of the order of the field frequency in small photonic lattices [20] and, up 
to our knowledge, it is possible to inscribe as many as sixty waveguides in fused silica 
with parameters related to Glauber- Fock photonic lattices [5l], which are of the order 
required by our nonlinear model, with a present technical limit set in the two hundred 
inscribed waveguides. On this account, caution must be exerted in finding the adequate 
balance between lattice size, parameters of the nonlinear Rabi model, and initial state 
to be propagated. 

One of us has discussed somewhere else [39] how a tight binding model of photonic 
lattices can be solved via matrix methods, idtS^"^^ = H^'^^E^'^^ leading to E^^^^t) = 
g-«//(±)t^(±)^g'j^ as long as the functions d^'^^j) and /(j) are time independent. In 
short, the dispersion relation of the truncated parity optical lattices of size N are given 
by the roots of the characteristic polynomial p^^\X) = described by the three-term 



recurrence relations: 

pit\X) =1, (18) 

#)(A) = [S^\0)-X]plt\X), (19) 

Pst^A) = [d^^\2k-l)-X]p^^l,-glf'i2k)p2k-2, (20) 

p£|i(A) = [d^^\2k) - X] pgl, - glf'{2k)p,,.u k > 1. (21) 



Note that in the nonlinear JC case, g^ = 0, its truncation depends heavily on the 
initial state of the field thanks to the fact that the Hamiltonian conserves the number 
of excitations. In the nonlinear Rabi case the truncation depends heavily on both the 
initial state and the value of the coupling parameter g+, it increases rapidly with the 
value of g+. For experimental realizations this just means that the lattice must be large 
enough to keep the propagated classical field far from the last segment of waveguides. 
Typical measurements of interest in the quantum optics community include, the 
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{n{t)) = {iP^+\t)\n\ij^^\t)) + {ij^-\t)\n\^^-\t)) 

|2 



Af-1 



6 

(22) 
(23) 



i=o 

which is nothing else than the sum of the output intensity at each waveguide weighted 
by the waveguide position; this quantity has been called the average center of mass of 
the propagating intensity in other context [5TH53] . Another relevant quantity is the 
mean atomic excitation energy, 

(aM) = {^^^\t)\a,\^M^t)) + {i;^~\t)\a,\^(~\t)), (24) 

(Af-l)/2 {N-l)/2 



\c{+) 
\^2j 



j=0 



\^2j+l 



j=0 



it)? 



(25) 



which is the addition of the intensity on the odd (even) waveguides minus the intensity 
on the even (odd) waveguides in the positive (negative) parity lattice. The atomic Berry 
phase is other quantity of interest, 

{a,{t)) = {m\^^\m), (26) 



Af-1 



(27) 



which clearly requires knowledge of the relative phases of the field at the outputs as it 
correlates the field amplitudes of positive and negative parity lattices. The fidelity of 
the time evolved state with respect to the initial state is another quantity of interest 
that requires knowledge of the relative phases of the amplitudes at the output of the 
parity lattices, 

^=|(^(0)|^(t))|, (28) 



Af-1 

$:(£r'(o))-£f'w + (f;-'(o))-f;-'w 

j=0 



(29) 



Finally, in order to recover information about the mixedness of the reduced atomic 
system, it is possible to calculate von Neumann entropy, 

S = -Tipa In Pa (30) 

where the reduced atomic density matrix is given by 

Pait) = T:Tf\m){mi 




(+) 12 
2i+l I 



c{+) c{-) 
\^2j 



c{~) (c(+) 
^2j+l l^2i+l 



c{+) lev- 
^27 + 1 \ ^2i 



'{-) 
'2j+l 



c(-) (c{+) 
^2j [<^2j 



l^2j I 



(+)|2 

2i I 



where the time dependence has been obviated for the sake of space, Sf" = £tit). 
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Now, the most common initial state configurations in quantum optics involve 
number or coherent states. We can simulate a number state times an atomic 
superposition, 

\^Pf{0)) = Ce\n, e) + Cg\n, g) , (33) 

Cg|+,n) + Ce|-,n), n even, 
Ce|+, + Cg|— , n), n odd. 

That is, just one input port is used in each parity lattices. In the case of n even the 
input at the positive and negative parity lattices are = \cg\ and si ^ = e'^\ce\ with 
the phase given by the phase difference between Ce and Cg and the equivalent for n 
odd. For coherent states, say the simplest case, 

\M0)) = \a,9), (35) 

(36) 

+, 2n) + —= 1-, 2n + 1) I . (37) 

V(2n + 1)!' " 



n=0 

This state is simply obtained by impinging a field on the 0th waveguide of a Glauber- 
Fock lattice [39|IM] and propagating it for the corresponding distance. Then, the output 
at the even waveguides is sent to the even waveguides of the positive parity lattice and 
that of the odd waveguides to the odd waveguides of the negative parity lattice. A 
highly entangled state involving coherent states can easily be constructed by sending 
the output from the Glauber- Fock lattice to just the positive or negative parity lattice; 
for example, if we consider just the positive lattice: 

°° n 

^-\a\/2j2^\+^n) = \a+,g) + \a_,e), (38) 




n=0 



where the even (odd) coherent state is defined by \a±) = \a) ± \ — a). The Glauber- 
Fock lattice can be used to simulate input corresponding to both coherent and displaced 
number states. 



4. An example: The Buck-Sukumar model 

In order to present a practical example, let us refer to the BS model, 

Hbs = ^ffi + + 9 (^SV^ a+ + y/h a^(T_ j , (39) 

that describes the interaction of a two-level system with a field under intensity dependent 
coupling. It is related to our model Hamiltonian (jj]) by setting = 0, g- = g, f{n) = n 
and h{h) = Ufh. Thus the positive and negative parity lattices, that classically simulate 
the BS model, are described by the coupled differential sets in Eq. (fT4l) - (fT7l) with the 
photon number functions /(n) = h and h{h) = ujfh and coupling parameters g_ = g 
and g+ = 0. 
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Figure 1. The classical simulation of the time evolution for the separable initial 
state IV'(O)) = IcK+iff) with a ~ 5 under BS dynamics on resonance, loq ~ Uf, and 
coupling parameters g_ = O.lojy, g+ = 0. (a) Propagation of the initial field in 
the corresponding positive parity photonic lattice of the classical simulator. The time 
evolution of the (b) mean photon number, (c) mean atomic excitation energy, (d) mean 
von Neumann entropy, and (e) fidelity reconstructed from the classical simulation. The 
lattice is composed by three hundred coupled photonic waveguides. 

In Figs. [1] and 121 we show numerical results for the classical simulation of the BS 
model on resonance and coupling g = O.luf] i.e. we use a photonic lattice described by 
the differential sets f ir^lfTTjl with parameter values Uq = Uf, g_ = O.luf, and = 0. The 
numerical propagation considers a photonic lattice composed by three hundred coupled 
waveguides. We consider the initial state |'?/'(0)) = \a^,g) with parameter values a = 5 
in Fig. [H This is a separable state with positive parity; i.e. just the positive parity 
photonic lattice is needed to classically simulate its evolution. The intensity of the 
light field is shown in Fig. [U^a). The time evolution of the mean value for the photon 
number, atomic excitation energy and von Neumann entropy are shown in Figs. [U^b)- 
[U^d). Figure [U^e) shows the time evolution of the fidelity, note how the evolution of this 
separable initial state returns periodically to its original state. In Fig. O we show the 
classical simulation of an entangled initial state \ip{0)) = g) + \a_,e) with identical 
parameter values as those in Fig. [1] Again, the state has positive parity and just 
the positive parity photonic lattice is needed to simulate the quantum system. The 
evolution of this entangled initial state also returns to its original state as witnessed by 
the fidelity in Fig. Et^e). 

For the sake of giving a complete example, we will also consider counter rotating 
terms, 

H = ujfh + + g- (a^fh (T+ + \fh <x^a\ + g^ (a\fh a_ + \fh a^a+ j 



A classical simulation of nonlinear JC and Rabi models 



50 
25.05 

(n) 

24.95 






(b) 



40 



t (units of ujf) 




t (units oi u!f) 



Figure 2. The classical simulation of the time evolution for the entangled initial 
state I'lpiO)) = \a+,g) + \a^, e) with a = 5 under BS dynamics on resonance, ujq — ujf, 
and coupling parameters g- = 0.1a;/, g+ = 0. (a) Propagation of the initial field 
in the corresponding positive parity photonic lattice of the classical simulator. The 
time evolution of the (b) mean photon number, (c) exponential of the mean atomic 
excitation energy, (d) mean von Neumann entropy, and (e) fidelity reconstructed from 
the classical simulation. The lattice is composed by three hundred coupled photonic 
waveguides. 



(40) 

We consider a model on resonance with parameter values ujq = Uf and g- = Qj^ = 2uf. 
Figure [3] shows the numerical results for the propagation of a light field that simulates 
the initial state \ip{0)) = |0,e) which has negative parity and corresponds to impinging 
the first waveguide of the negative parity photonic lattice. Again, the intensity of the 
light field is shown in Fig. [3t^a). The time evolution of the mean value for the photon 
number, atomic excitation energy and von Neumann entropy are shown in Figs. EJ^b)- 
|3|^d). The numerical propagation considers a photonic lattice of size two thousand and 
the probability of finding light at the last waveguide has a maximum value of 7 x 10~^ 
within the parameter range considered here. Note that this is a thought experiment at 
the time because current technology can produce a couple hundred coupled waveguides 
at most. 

5. Conclusion 

We showed how the parity of the nonlinear Rabi model in ([2]), which is unphysical for 
most parameter sets with current quantum technology, can be exploited to simulate its 
dynamics by propagating classical light fields in a set of two photonic lattices, each 
one corresponding to one of the parity subspaces. We discussed how initial states 
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Figure 3. The classical simulation of the time evolution for the separable initial state 
|'(/;(0)) = |0, e) under BS plus counter rotating terms dynamics on resonance, = '^f^ 
and coupling parameters g- = g+ = 2a;/. (a) Propagation of the initial field in 
the corresponding negative parity photonic lattice of the classical simulator. The time 
evolution of the (b) mean photon number, (c) mean atomic excitation energy, (d) mean 
von Neumann entropy, and (e) fidelity reconstructed from the classical simulation. The 
lattice is composed by two thousand coupled photonic waveguides. 



in the quantum system involving Fock states map to impinging the photonic lattices 
at specific waveguides. Initial states involving coherent or number displaced states 
can be engineered by propagation through Glauber-Fock lattices. As an example, we 
presented numerical results for the classical simulation of the Buck-Sukumar model 
with initial states within the positive parity subspace describing fully separable and 
maximally entangled states. Also, we considered the Buck-Sukumar model including 
counter-rotating terms for an initial state in the negative parity subspace. 
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